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Abstract 

In this paper, we consider multiple channels and wireless nodes with multiple transceivers. Each node 
assigns one transmitter at each available channel. For each assigned transmitter the node decides the 
power level and data rate of transmission in a distributed fashion, such that certain Quality of Service 
(QoS) demands for the wireless node are satisfied. More specifically, we investigate the case in which the 
average SINR over all channels for each communication pair is kept above a certain threshold. A joint 
distributed power and rate control algorithm for each transmitter is proposed that dynamically adjusts 
the data rate to meet a target SINR at each channel, and to update the power levels allowing for variable 
desired SINRs. The algorithm is fully distributed and requires only local interference measurements. 
The performance of the proposed algorithm is shown through illustrative examples. 

I. Introduction 

Throughput is an important metric in wireless networks and its improvement is achievable by 
efficient use of the wireless channel. Since each transmission corresponds to a spatiotemporal 
propagation of radio waves that are received by all nodes in proximity utilising the same channel, 
nodes interfere with each other when they use the same channel simultaneously. This is called co- 
channel interference (CCI) |[T|. In addition, the power of each transmitter in a wireless network 
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is directly related to the resource usage of the link and it is a valuable resource, since the 
batteries of the wireless nodes have limited lifetime. Increased power ensures longer transmission 
distance and higher data transfer rate. However, power minimisation not only increases battery 
lifetime, but also increases the effective interference mitigation which in turn, increases the 
overall network capacity by allowing higher frequency reuse. Furthermore, the near-far problem 
Q is commonly solved by dynamic power adjustments by the transmitters. Dynamic power 
control in wireless networks allows devices to setup and maintain wireless links with minimum 
power while satisfying certain constraints on QoS. As a result, power control (also referred to as 
transmission power control) has been a prominent research area with increased interest during 
the last two decades. 

In this paper, we combine power control with transmission rate in the presence of multiple 
channels. The main goal for each communication pair is to retain the average of the SINRs 
(or the transfer rates) from all the available channels above or equal to a certain threshold. A 
distributed scheme is developed in which individual users can cooperate in such a way that 
otherwise infeasible states can be achieved in a distributed manner without the need of interlink 
communication among users and centralized computation as required in a centrally controlled 
wireless environments. In order to be able to reach states that are infeasible to be achieved by 
simultaneous transmissions, it is necessary to allow for variable desired SINRs to be assigned 
in each of the available channels. Hence, we relax the hard constraints for SINRs and try to 
achieve the desired QoS on average by making use of multiple channels. As a result, the problem 
targeted in this paper is how communication pairs can achieve on average the QoS targets that 
belong in the convex hull of the system's feasible set, but would probably be impossible to be 
achieved with the network setting by simultaneous transmissions. The difference of the problem 
being solved in this study compared to other related studies (e.g. for opportunistic transmission 
and Multi-Input Multi-Output (MIMO) interference systems) relies on the fact that we consider 
different frequency channels in which the network configuration might be very different, and also 
consider the general SINR regime in which the problem cannot be transformed into a tractable 
convex optimisation problem. As a result, even though a considerable amount of work has been 
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done to characterize the capacity of the system and to find the maximum throughput of the 
system when multiple channels are available, to the best of our knowledge there is no work on 
the joint power control and transmission rate allocation for multiple channels, when multiple 
users appear in the network and they need to fulfil a QoS requirement in the general SINR 
regime. 

Early work in the field of power control for wireless networks [|2|, [[3| proposed power 
balancing, which equalizes the Signal-to-Interference Ratio (SIR) in all the wireless links. 
These algorithms need global information about the network setting. This capacity improvement 
initiated extensive research on power control with focus on the design of distributed algorithms 
(working synchronously [Q, ||5|, asynchronously [[6), [|7|, under constraints [8], with active link 
protection Q, in the presence of time delays [[T0|-|[T4| etc.) to meet a prefixed SINR target (hard 



constraint), determined by the QoS requirements. The prefixed SINR target tracking (which is 
the condition for inelastic traffic) ensures that a constant transmission rate can be sustained. If 
a feasible solution exists, then there exists a unique solution that minimizes transmit power in a 
pareto sense. But, if not, then the performance of the whole network degrades and the capacity 
is deteriorated. The target tracking approach is suitable for real-time, delay-sensitive applications 
like mobile phone services. 

In view of the proliferation of wireless data though, it is essential to investigate further 
transmission schemes, that is, techniques that facilitate elastic and/or opportunistic traffic should 
be considered, where time-varying rates are allowed and large delays are tolerated, such as in 
the World Wide Web (WWW) and video streaming. When each node is assigned its desired 
SINR, it has no knowledge of the network and as a result, the combination of all users' QoS 
requirements are sometimes impossible to be fulfilled in a single channel. If nodes could exchange 
information and obtain full information about the network setting, then they could adjust their 
desired SINRs so that the network becomes globally asymptotically stable when all transmitters 
operate simultaneously. Since every node affects all the other nodes in the network, it is not 
possible (especially in big networks) to acquire the knowledge required about the whole network. 
Furthermore, even in small networks consisting of two communication pairs only (i.e. four 
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nodes), simultaneous transmissions for the desired SINRs may not be feasible. 

In general, the fluctuation of wireless channels can be exploited using power and transmission 
rate control in order to meet QoS requirements, i.e., a node can increase its transmit power 
whenever the interference at its receiver is low and decrease it when the interference is high. The 
adaptation of transmission rate to the actual propagation and interference conditions is facilitated 
by modem adaptive modulation and coding schemes. The transmission rate can be selected so 
that the Bit Error Rate (BER) is sufficiently small and this can be achieved by adjusting the 
transmission rates according to the SINR. In that way, more information is transmitted when the 
channel conditions are favourable by adjusting the transmission rate accordingly. This approach 
enables the improvement of the network's convergence and the satisfaction of heterogeneous 
service requirements. 

It is well known that the adaptation of the transmit power, data rate, and coding scheme 
increases spectral efficiency. The IEEE 802.11b scheme allows nodes to increase their transfer 
rate up to 11 Mbps, depending on the SINR at the receiver. The performance achieved through 
power control can be further improved by allowing for dynamic adjustment of the transfer 
rate based on the SINR. There exist some approaches in the literature that deal with dynamic 
adjustment of the desired SINR, and hence the data rate. Some consider a single channel only 



(e.g. [ 15 1-[ 18 1) whereas others consider the multi-channel case (e.g. [ 19 1-[21 1) as well. Concepts 
and related work on both approaches are briefly described in turn. 

Firstly, we describe the approach that supports the use of adaptive desired SINR for the single 



channel only. Some schemes (e.g. p7| , p2| and p8| ) tried to deal with the case in which 
the desired SINR could be adaptive, depending on the channel conditions, but the main idea is 
restricted within the limits of a compromise on the QoS demand. Their approach enables the 
improvement of the system convergence and the satisfaction of heterogeneous service require- 



ments. In [23 1, [24 1 an opportunistic power control algorithm is proposed which stems from the 
concepts of opportunistic communications [ [25| , [ [26| . This algorithm provides tunable parameters 
in order to exploit the trade-off between throughput and power consumption. The algorithm is 
proven to converge to a unique fixed point. However, it does not consider a QoS requirement. 
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but rather tries to maximize the throughput for individual users. The same problem has also 
been targeted by p7| , in which a class of distributed power control algorithms that exploit the 
channel variability opportunistically is presented. These algorithms are restricted to maintain the 
desired QoS in a certain range, which is specified by the node. Hence, these algorithms can be 
converted to a conventional algorithm that requires a prefixed QoS requirement. Nevertheless, 
none of the approaches associate the QoS requirement with a cost or utility function that leads 
to a solvable power control problem. Consequently, there is no way to guarantee that on average 
their QoS targets are fulfilled. In other words, there is no way to maintain a constant average 
data rate. 

Since wireless nodes may use different channels and they have different locations at different 
times, resource allocation can take advantage of the diversity in space, time and frequency. 
Wireless technology standards provide a radio-frequency (RF) spectrum with a set of many 
non- overlapping channels and a node has the option to choose on which channels to transmit. 
Therefore, we can make use of extra channels, so that nodes experiencing much interference by 
other nodes can choose to operate in a different channel, if possible. However, since it is possible 
to establish diversity by introducing multiple transceivers or transceivers with multiple antennae 
on a single node, communication between two wireless nodes should not be confined on a single 
channel only, but make use of all the available channels. Multi-Input Multi- Output (MIMO) links 
are antenna arrays at both ends of a link and transmit parallel streams at the time and frequency 
channel. Channels arising from the use of spatial diversity at both the transmitter and the receiver 
have been considered in the literature (e.g. [28] and reference there in) and it is shown to have 
improved capacity with respect to Single-Input Single-Output (SISO) systems because of the 
use of parallel channels; such a mutli-access network with MIMO links has been considered 



in previous studies (such as [29|-[31| and references therein). As aforementioned there is a 
significant difference of the problem being solved in this study and MIMO interference systems. 
In this study, we consider different frequency channels in which the network configuration might 
be very different and also consider the general SINR regime in which the problem cannot be 
transformed into a tractable convex optimisation problem. Even though a considerable amount 
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of work has been done to characterize the capacity of the system and to find the maximum 
throughput of the system when multiple channels are available, to the best of our knowledge 
there is no work on the joint power control and transmission rate allocation for multiple channels 
when multiple wireless nodes appear in the network and they need to fulfil a QoS requirement 
in the general SINR regime. 

The rest of the paper is organized as follows. Section|n]establishes the notation used throughout 



the paper. Section III describes the mathematical system model. Section IV presents some 
preliminary results on the conditions for feasible networks and provides sufficient conditions for 
stability of the Foschini-Miljanic (FM) algorithm [|5| using Lyapunov Stability theory, a useful 



preliminary result for the subsequent analysis. Section V-A addresses the feasibility regions based 



on the channel model considered. Section V-B formulates the problem being targeted. Section 



V-C derives the main result of this work, in which we develop a distributed scheme that allows 



variable desired SINRs and facilitate the use of multiple channels. Section VI demonstrates the 



validity of the proposed scheme through illustrative examples. Finally, Section VII concludes 



the paper with a brief discussion on the algorithm developed as well as future directions. 

II. Notation 

The sets of complex, real and natural numbers are denoted by C, IR and IN, respectively; their 
positive orthant is denoted by the subscript + (e.g. C+). Vectors are denoted by bold letters 
whereas matrices are denoted by capital letters. and A^^ denote the transpose and inverse 
of matrix A respectively. For two symmetric matrices A and B, A >- {'^)B means that A — B h 
(semi-)positive definite. By / we denote the identity of a squared matrix. \A\h the elementwise 
absolute value of the matrix (i.e. \A\ = [\Aij\]), A{<) < B is the (strict) element-wise inequality 
between matrices A and B. A matrix whose elements are nonnegative, called nonnegative matrix, 
is denoted by A > and a matrix whose elements are positive, called positive matrix, is denoted 
hy A > 0. a{A) denotes the spectrum of matrix A, \{A) denotes an eigenvalue of matrix A, 
and p{A) denotes its spectral radius. det{A) denotes the determinant of a squared matrix A and 
diag(xi) the matrix with elements xi, X2 , . . . on the leading diagonal and zeros elsewhere. 
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Further notation used in the paper is tabulated below (Table [T]): 



Notation 1 Notation used in the paper: 



M 


The 


C 


The 


g 


The 


r 


The 


n 


The 


c 


The 


9ij 


The 


Pi,k 


The 


Wi,k 


The 


■^i,k 


The 


^i,k 


The 


li 


The 



:th 



in. System Model 

The system model can be divided into two levels: level 1 describing the network as a 
whole; and level 2 describing the channels. Thus, we have the network model and the channel 
model. At the network level, the model concerns the general topology of the nodes and their 
characteristics. At the channel level, the model describes the assessment of the link quality 
between communication pairs and the interaction between the nodes in the network. 

A. Network Model 

Consider a network where T denotes the set of transmitters and IZ denotes the set of receivers 
in the network. The links are assumed to be unidirectional and each node is supported by omni- 
directional antennae. At each time instant, each node can act as a receiver or a transmitter only 
due to the half-duplex nature of the wireless transceiver. Each transmitter aims at communicating 
with a single node (receiver) only. 

B. Channel model 

As mentioned above, the link quality is measured by the SINR. The channel gain on the link 
between transmitter i and receiver j is denoted by gij and incorporates the mean path-loss as 
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a function of distance, shadowing and fading, as well as cross-correlations between signature 
sequences. All the gij's are positive and can take values in the range (0, 1] (see Figure [T|). The 
power level chosen by transmitter i is denoted by pi and the intended receiver is also indexed 
by i. Ui denotes the variance of thermal noise at the receiver z, which is assumed to be additive 
Gaussian noise. 

The interference power at the i\h receiver consists of both the interference caused by other 
transmitters in the network ^^^7- QjiPj (where Tli denotes all the transmitters j in the network 
that interfere with transmitter's i communications, i.e. j i,j E T), and the thermal noise z/j 
in node z's receiver. Therefore, the interference at the receiver i is given by 

= ^ gjiPj + Vi, (1) 



while the SINR at the receiver i is given by 



(2) 



Eier-, 9jiPj + ' 

Due to the unreliability of the wireless links, it is necessary to ensure QoS in terms of SINR in 
wireless networks. Hence, independently of nodal distribution and traffic pattern, a transmission 
from transmitter i to its corresponding receiver is successful (error free) if the SINR of the 
receiver is greater or equal to 7^ (Fj > 7j), called the capture ratio which depends on the 
modulation and coding characteristics of the radio. Therefore, 

> 7i (3) 



Remark 1. We have not specified any model for determining the positions of the nodes, since 
we investigate the general case of a network that any position could be possible. Further, we 
have also not specified any model related to the propagation of signals. In our context, these 
two models will ultimately specify the channel gains, Qij. Nevertheless, they are of secondary 
importance in this study since it is focused on how the QoS is improved given the channel gains 
and it only depends on the power of the received signals. Note also that the effect of nodes' 
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Fig. 1. An example of a network consisting of two communication pairs only. Each pair i consists of a transmitter Si and 
a receiver Ri connected with a solid line while the grey dotted arrows indicate the interference that transmitters cause to the 
neighboring receivers. 



mobility is not consider in this study. However, this could be relaxed to the case of low mobility, 
where the link structure is expected to change slowly with respect to the packet rates and network 
updates. 

IV. Preliminaries 

A. System Feasibility 

Inequality ([3]) which depicts the QoS requirement of a communication pair i while transmission 
takes place is equivalent to the following condition: 

In matrix form, for a network consisting of n communication pairs, this can be written as 

P > TGp + 77, (5) 
where r = diag(7i), p = r ... p^ \ , i]i = 'jiUi / gu and 




, if « = J, 
, if ^ ^ j- 
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Let C = TG such that 



then (|5]) can be written as 



, if « = J, 
= \ (6) 



(/ - C)p > rj. (7) 



The matrix C has strictly positive elements off diagonal and it is reasonable to assume that is 
irreducible [ [32| , since we are not considering totally isolated groups of links that do not interact 
with each other. By the Perron-Frobenius theorem [j32|, we have that the spectral radius of the 
matrix C is a simple eigenvalue, while the corresponding eigenvector is positive componentwise. 
The necessary and sufficient condition for the existence of a nonnegative solution to inequality Q 
for every positive vector rj is that (/ — C)^^ exists and is nonnegative. However, {I — C)^^ > 



if and only if p{C) < 1 |33| (Theorem 2.5.3), [34|. Therefore, the necessary and sufficient 
condition for ^ to have a positive solution p* for a positive vector rj is that the Perron- 
Frobenius eigenvalue of the matrix C is less than 1. Hence, a network described by matrix C 
is feasible if and only if p{C) < 1. 



B. Stability of the FM algorithm using Lyapunov Theory 

Given a network, we know that a feasible solution to the power control problem exists if the 
Perron-Frobenius eigenvalue of matrix C is less than one (p(C) < 1), where matrix C is given 
by ([6]). If the condition holds, the FM algorithm [5] is asymptotically stable. Since the FM 
algorithm is linear, when the stability condition holds we know that there exists a quadratic 
Lyapunov function which establishes the stability of the system. 

Now, the notions of Lyapunov and D-stability are being used for the derivation of some prelim- 
inary results on the stability of the FM algorithm, defined by the following differential equation 
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0: 

where ki E R+ denotes the proportionality constant and 7^ denotes the desired SINR. It is 
assumed that each node i has only knowledge of the interference at its own receiver. 
By defining e(t) = p(t) — p*, where p* denotes the vector of power levels when the system is at 
the equilibrium, the stability of the system will be assessed. With the above system of equations 
in mind, we note the following important result that is central to the conclusion of this paper. 

Lemma 1. Suppose the spectral radius of matrix C in ([6]) z^' less than 1, then the quadratic 
Lyapunov function 

V{e) = e^A^DAe, (9) 
where A = KH, K = diag(A;i), H = I — C and D is a positive diagonal matrix, stabilizes ([8]) 

for ji, Qji, z/j > 0, for any initial state Piif)) > and for any proportionality constant, > 0. 



Proof: The power control algorithm (10) can be written in matrix form as 



p{t) = ~KHp{t) + Kri (11) 



where K = diag{ki) and 




H^i = < (12) 



Let A = KH, q = Kr] and e(t) = p(t) — p* ^ e{t) = p{t). Then, ( [TT] ) can be expressed as 



p(t) = -Ap(t) + q. (13) 
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At equilibrium, p = and hence Ap* = q. Therefore, ( [13] ) is written as 

Pit) = e{t) = -Apit) + Ap* = -A(p(t) - p*) = -Ae(t), (14) 

giving e(t) = —Ae{t). Consider the candidate Lyapunov function 

V{e) = e^A^DAe, (15) 

where D is a positive diagonal matrix. V'(e) > V e(t) 7^ and F(e) = for e = 0. In addition, 
F(e) is radially unbounded, i.e., V{e) — )■ 00, when ||e|| — )■ 00. Therefore, if V{e{t)) < 0, V 
e(t) 7^ 0, then the equilibrium point is globally asymptotically stable. 

V{e{t)) = e^{t)A^DAe{t) + {t) A^ D Ae{t) 

= {-Ae{t)YA^DAe{t) + {t)A^ DA{-Ae{t)) 
= ~e^{t)A^{A'^D + DA)Ae{t) 

From the equation above V{e(t)) < 0, V e(t) 7^ 0, if and only if A^D + DA y 0. Since A is a 
nonsingular M-matrix, there exists a positive diagonal matrix P such that the matrix A^P + PA 



is positive definite [35| (condition H24,, p. 136). As a result, there exists a positive diagonal matrix 
D for which A^D + DA >- and therefore ( [T5| ) is a Lyapunov function that stabilizes the system 



given by (11 1. 



V. Joint Distributed Power and Rate Control 
A. System Feasibility Region 

In this subsection, we will provide some intuition on the SINR (and capacity) regions of a 



network, based on the model described in Section III More specifically, we show that they are 
non-concave sets with respect to the transmitters' power and we suggest that it is possible to 
reach regions within the convex hull of these non-concave sets. Note that previous works, such 



as [36 1, converted the capacity regions to concave sets by working in the high SINR regime. 



and hence assuming that the capacity is given by log2(rj), instead of log2(l + Fj), where Fj is 
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given by formula ([2]). 



Definition 1. (Feasibility) /|37|/. A set of target SINRs Fj is said to be feasible with respect 
to a network, if it is possible to assign transmitter powers pi > so that the requirement in 
inequality ([3]) is met for all nodes transmitting in the network. Analogously, the power control 
problem is said to be feasible under the same conditions. Otherwise, the target SINRs and the 
power control problem are said to be infeasible. 

We consider a wireless ad hoc network with M transmitters operating in a single channel. Let p 
be the M-dimensional power vector with pi the power of transmitter i in the network. Moreover, 
we include upper bounds on the power level, since nodes operate with small batteries which 
impose limitations on the transmission power. Let Pmax be an M-dimensional vector that contains 
the maximum power level for each transmitter in the network. Note that pi lies in the closed 
interval [0,Pi,max\- We define 11 to be the set of all possible powers levels in the network, namely, 
n(p) — {p : < p < Pmax}- By definition, 11 is described geometrically by a rectangle for 
two transmitters only and an M-orthotope for M transmitters. Therefore, the set of powers is 
always a convex set. 

We denote by Fj the SINR of transmitter i given by (|2]), and Ti ^ax its maximum SINR, defined 
by 

^i,max = {^i ■ Pi = Pi,max , = Vj ET,] ^ i} . (16) 

We denote by T the Af -dimensional vector with Fj the SINR of transmitter i. We define the set 
to be all feasible data rates that can be achieved in the network at any instant by simultaneous 
transmissions, i.e., ^(F) = {F : p G H}. Since Fj's are not concave functions of the transmitters' 
power levels, depending on the network geometry and topology, ^ is not necessarily a convex 
set. Consequently, we define the set $ to be the convex hull of ^. i.e., $ = Co(^) . 
According to the definition of Fj and the assumption that there exists thermal noise in the network, 
the set $ is always closed and bounded. Hence, although the SINR formula is generally a non 
convex function, the convex hull operator transforms the SINR region, $, into a convex set. 
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Formally, based on the above definitions, we define the SINR region of a wireless network as 
the convex hull of all the feasible SINRs. 

If $ — \E' = 0, then \1' is a convex set, which means that the system performance cannot be 
improved further. If, for a chosen T there exists a solution to the system for simultaneous 
transmissions, then F G \1'. However, there exist cases where F G $ — ^. In such cases, a 
scheme that partitions time or introduces multiple channels should be introduced. 
By Caratheodory's theorem [38|, if there is a point x G R^^, lying in the convex hull of the set 
\E', there exists a subset ^' of ^ consisting of no more than M + 1 points such that x lies in the 
convex hull of Therefore, any point F G R*^, that lies in $, can be achieved with a linear 
combination of no more than M + 1 points in ^. This suggests that, by efficiently partitioning the 
time between different states, target SINRs can be achieved that would otherwise be infeasible. 
Also, the solution consists of at most 5'+ 1 different states, where S is the number of transmitters 
in the network. 

Remark 2. The same arguments can be used for the case where the capacity Rj is considered, 
due to the monotonic relation with SINR (since Rj = log2 (1 + Fj) - Shannon's capacity). Again, 
capacities are not concave functions of the transmitters' power levels. They depend on the 
network geometry and topology, and hence, the capacity region of a wireless network is defined 
as the convex hull of all the feasible capacities in the network. 

Some representative examples of networks are presented in Figure [2] in order to show the ideas 
aforementioned. Particularly, we show that there exist networks for which it is possible to attain 
data rates that would be otherwise impossible, with time division or simultaneous transmissions 
over many channels. In the following examples (Figure [2]), on the left-hand side the wireless 
ad-hoc network is depicted consisting of the communication pairs. Each of these pairs consists 
of a transmitter and a receiver. The transmitters appear as nodes with an S and the receivers 
appear as nodes with an R. Communication pairs have the same index and are connected by a 
solid line {Si — > Ri}. The grey dotted arrows denote the interference invoked and it is shown to 
some of the connections for illustration purposes only. On the right-hand side, the feasible data 
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(a) A wireless ad-hoc network is depicted, consisting of 
n = 4 nodes and hence two communication pairs. 




(b) In region K, simultaneous transmissions can take place 
(C G 4*), whereas in region L, a higher throughput can be 
achieved by partitioning the time between different states of 
data rates (C G "l? — '5)). In region M, there is no feasible 
throughput for the network considered (C ^ 




\ 912 



\f/22 



(c) A wireless ad-hoc network is depicted, consisting also 
of n = 4 nodes and hence two communication pairs. 
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(d) The sets of feasible data rates (4' C "I") for this network 
are illustrated. In order to achieve the set of data rates in 
region L, one of the transmitters switches to a different 
mode of switching between states, whereas the other keeps 
transmitting at its maximum power. 




(e) A wireless ad-hoc network is depicted, consisting of 
n — 6 nodes and hence two communication pairs. 



(f) The two pairs for which the receiver-transmitter distance 
is small can transmit simultaneously and $ — = 0, 
whereas the third, due to the large distance between the 
transmitter and receiver has a low throughput and it deteri- 
orates when any of the other pairs transmits simultaneously. 

ptgolji- 3t)^Qhfe left-hand side the wireless ad-hoc network is depicted consisting of the communication pairs. The transBBtyfS 
appear as nodes with an S and the receivers appear as nodes with an R. Communication pairs have the same index and are 
connected by a solid line {Si — > Ri}. The grey dotted arrows denote the interference invoked. On the right-hand side, the 
feasible data rates of the communication pairs along with the capacity region of the network are illustrated. 
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rates of the communication pairs along with the capacity region of the network are illustrated. 

B. Problem formulation 

We now formulate the problem of distributed power control in the case which we have more 
than one channels available and the average SINR of the system is required to meet certain QoS 
specifications. Therefore, we have a variable SINR threshold in each channel, while the average 
of all the SENR targets should meet the overall required QoS. 

Let i e T in channel k E C, where T and C are the sets of wireless transmitters and available 
channels, respectively. As before, 7, denotes the desired average SINR of user i for all 
channels. We introduce Xi^k to be the allotted (or desired) SINR target for user i and Wi^k is the 
effective interference at receiver i in channel k, and is given by 



Therefore, for each channel k we require that the SINR is greater than or equal to Xi^k- 

— > Xi^k, Xi^k > (18) 

Wi,k 

where Pi^k is the power of user i in channel k The target SINR (7^) should be met on average 
over all N channels, i.e.. 



fceC 

Therefore, we are required to develop an algorithm that adapts the targeted SINR in each channel, 
accounting for the individual channel's conditions, in addition to updating the power. 

C. Main results 

We propose a distributed algorithm that not only updates the power level ^ of each node i in 
a certain channel k, but also updates the desired SINR Xi^k for each channel in order to take 
advantage of the channel conditions. We show that this algorithm is able to find feasible solutions 
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that would be inaccessible via conventional power control algorithms (e.g., the FM algorithm) 
that operate in a single channel only. The proposed algorithm aims to find an equilibrium for 
the system consisting of all the communication pairs making use of all the available channels, 
if a feasible solution exists. 

Definition 2. A feasible solution exists if there exist SINK targets = [ xi k X2k ■ ■ ■ xn k V 
such that 

p(C(x,))<l, \/keC (20) 

and 

kec 



(21) 



In the following theorem we state that the proposed distributed algorithm is asymptotically stable 
and converges as long as the allotted SINRs belong to a feasible solution. 

Theorem 1. Assume that 7^, gji, Vi G R+, for any initial state Pi(0) > and for any propor- 
tionality constants Ci^k,bi,k,Ci ^ '^+' <^^d. suppose the spectral radius of matrix C^'K^it)) V t in 
(|6]) remains less than 1. Then, the distributed power control formula 

and the distributed allotted SINK target update formula 

Xi,k{t) = hi^t) ^ ( 7i - ^ X] ^i^^) ) ~ '^i,kit) iXi,kit)WiAt) - Pi,kit)) (23) 

\ fcec / 

converge to a feasible solution. 

Proof: For each transmitter i E T, WQ let the utility function be 

Ui := did I 7i - ^ XI ) + X] t*^*'^ {xi,kWi^k - Pi,k)f (24) 
\ kec / kec 
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Pm ] ^ '^k — \ x\ X2 ... X 



M 



, D 4 



where q fc, di and Q are positive constants. 

Define efc = - p^, ^ (yp\ v\ . 
diag((ii) and ^(x^) is the same as before, but with Xi^k in place of 7^ for each node respectively. 
Explicitly, since A=KH, we express H as 



Hi 



,ifi^ j, 



(25) 



-Xi,kf: , if ^ 7^ J- 



Therefore the utility function for the whole network, with M communication pairs, is given 



by, 



M 



did ( 7j - ^ X] Xi^k I [ci,k {xi,kWi,k - Pi,k 

\ k&C / keC 

M / ^ \ 2 M 

i=l \ keC / i=l keC 

M / ^ \ 2 M 

i=i \ fceC / fcec 1=1 

i=i \ kec J kec 



k {xi^kWi^k -Pi,k)y 



Xi,kWi,k - Pi,k) 



U2 



By definition, the utility function is nonnegative, i.e., 

U>0, Vejfc^O and ^^Xi,fc^7i Vi. 

Also, 



(26) 



[/ = for efc = and ^ Xi^k = li Vi. 



(27) 



feec 



In addition, [/ is radially unbounded, i.e., U ^ 00, when ||efc|| — > 00 or/and ||xfc|| ^ 00. 

In what follows, we find the conditions for which tf < 0, V 7^ and X^feec Xi,k/N 7^ 7^. That 
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IS, 



M 



dU v-^ dU dxi h v-^ dU deu 

— => > ^ + \ ^ < 0. 

dt ^ ^ dxi k dt ^ dek dt 

i=i kec ' kec 



(28) 



for all efc 7^ and Ylk^hk/N ^ 7j. In Lemma [T| we showed that the Lyapunov function 
(|9]) that proves stability for the FM algorithm is of the same form as f/|; namely, = 
elA^{xk)DA{xk)ek. Convergence in case the matrix A is time varying, A(xk{t)), is proven 



in [39|, [40|. Therefore, there exists positive diagonal matrix D such that A^D + DA > —Q 
for any positive definite matrix Q, and given that p(C(xfc(i(:))) < 1, such that. 



EdU dek ^ ^ . ^ 
— — < 0> V efc ^ 0. 



fcec 



dek dt 



(29) 



Ui is independent of e.k{t) and therefore f/ < 0, V efc(t) 7^ 0. Thus, the power update formula 



( [22j ) stabilizes the system. Now, we show that the update formula for Xj ^ establishes that 

M 

^ ' N 



EE 



i=l keC 



dxi k dt 



(30) 



kec 



For a single transmitter. 



dx 



i,k 



\ kec / 



+ 2diWi^k {xi,kWi^k - Pi,k) 



(31) 



kec 



kec 



-2d,J2 

kec 

-2di J2 



-E 



kec 



Xi,k - Wi^k [Xi,kWi^k - Pi,k 



Xi k 



kec 



N 



li 



kec 



Xi,k (32) 
(33) 
<0, (34) 



where ( [34] ) is obtained by substituting ( [23] ) into ([33]). Since Xj ^ does not appear in the other 
nodes' utility functions, it is shown that for each user, the allotted SINR target's differential 
equation makes the rate of change of the utility function of node i with respect to Xi^k negative. 
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We have proved the update formulae converge to a feasible solution under the condition p(C(xfc)) < 
1. Hence, the suggested algorithm converges to a feasible solution as long as the desired SINRs 
belong to the feasible set of SINRs. 

However, this might not be always the case, and p(C(xfc)) > 1 for some time-interval. In what 
follows, we introduce a condition on the update formulae such that the algorithm converges to 
a feasible solution, if one exists, even if p(C(xk)) > 1 over bounded time intervals. 



Theorem 2. Define 



Oiit) := 7i - J^^^iAt) 



where N is the number of available channels, Xi^k is the allotted SINR target for transmitter 
i E T in channel k E C and '~fi is the average capture ratio desired at receiver i. The distributed 
power control formula 



Pi At) 







Cjjfc I Pi,k 



' if 7^ and Xi^k = 0, VA; G C, 
, otherwise. 

(35) 



and the distributed allotted SINR target update formula 



N 



It 



kec 



(36) 



converge to a feasible solution, if one exists, for all positive constants Ci^k G Q G R+, and 
provided the proportionality gain hi ^it) G R+ of the allotted SINR target update formula (|36]) 
is appropriately chosen such that the inequality 

/ \ -1 2 



biA^) 



N 



k [Xi,kWi,k -Pi,k) 



kec 



> '^MmA^) 



(37) 
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is satisfied by all nodes i in at least one of the channels k ^ C and 



i,k 



< 0. 



(38) 



Proof: For each transmitter i, we let the utility function ( [24| ) as before. The global utility 
function is given by the summation of all the individual utility functions and can be written as 



M 



M 



U = Y,U^ = Y,d^ 



i=l 



1=1 



\ kec / k 



Let (j)i^k = Pi,k, then (f)i^k = Pi,k- Hence, 

M M 



i=l 



i=l 



^ dUi 



dxi. 



— — Vi.k 



.k&C 



H,k 



By substituting Ui as given in (39) into (40), 



1=1 fcec ^ ' ' 



This implies that f/ < if 



X: 



i,k 



+ 20i,fc0- < 0, V A; G C and V i G r. 



Therefore, by substituting ( [36] ) into ( [42| ), the following condition is derived: 

n 2 



h,k{t) 



N 



I 7i - ^ X] Xi,k j - Wi^k {xi^kWi^k - Pi,k) 



> 2p,^k{t)pik{t). 



(39) 



(40) 



(41) 



(42) 



(43) 



For ij fc 7^ 0, if Pi^k{t) < 0, then since pffc(t) > 0, inequality ( [43] ) holds for any proportionality 
gain bi^kit) > 0, whereas if Pi^kit) > 0, then h^kit) can be adapted so that ( |43] ) holds. 



For Xj A; = and if Pi,k{t) > 0, then the inequality is not fulfilled in channel k. However, as 
long there exists at least a channel, say m, which satisfies Xj „i 7^ 0, then the proportionality 
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gain hi^rn{t) can be adjusted such that 



hi 



i,k 



< 0. 



(44) 



while fulfilling inequality (43) as well 



For Xi^k = \fk E C and 6i ^ 0, then Pi^kit) = 0. Hence ( |36l ) becomes 



Xi,k{t) = bi,k(t) 



fcec 



and 7^ 0. 



(45) 



For 9i = 0, then ^ 7^ V/c G C unless pi^k{t) = for which equilibrium is reached. Hence, 
the condition on the power update formula establishes that f/ < in all cases, apart from the 
case in which a feasible solution is reached. ■ 

Remark 3. If the maximum power pj.max is reached, the algorithm continuous to update both 
formulae as before, and if Pi^kif) < 0, then the power is updated. Hence, the upper bound of 
the power does not affect the analysis of this work. 

VI. Illustrative Examples 

To demonstrate the essence of the results derived, we present some simple examples. We 
consider a simple wireless ad-hoc network consisting of n = A nodes (i.e. two communication 
pairs Si — > Ri). The network structure is shown in Figure [3} The network is symmetric, so initial 
powers and/or desired SINRs should be different for the pairs in order to observe differentiation 
in their actions. 



911 



Distance 



922 



Fig. 3. A wireless ad-hoc network of n — 4 nodes, consisting of two communication pairs {5*^ -> Ri}. 



By varying the distance between the two communication pairs, the SINR region of the network 
is changed. Thus, we can either choose a convex SINR region to study or a non-convex one. 
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Based on this network prototype we examine the performance of our algorithm in different cases 
of network configuration. 

The parameters used for the cases studied are summarized below: 



Parameter 


Value 


Average SINR target (ji) 


3 


Constant (Q) 


20 


Proportionality constant (ki) 


1 


Initial Proportionality gain (6i_fc) 


200 


Noise 


0.04 fiW 


Maximum Power (p™,?^) 


10^ fiW 



TABLE I 

Parameters of the algorithm and the wireless networks used in the following examples. Power and 

NOISE are measured IN Watts (W) AND DATA RATE IN BITS PER SECOND (bits/s). 



The initial power levels (pi,k) for each node i and for each channel k are chosen at random. 
The upper bound for the maximum power for the simulations is set to lO^mW, large enough to 
allow the distributed algorithms to operate without this extra condition on the power, since we 
have not considered an upper bound in our derivations. The proportional gain 6j ^ is initially set 
to 200 for all transmitters in the network, which is high enough to guarantee that condition ( |37| ) 
is fulfilled and it is increased whenever required, depending on the condition so that it maintains 
a safe margin from the minimum value. 

A. Example 1: Feasible network 

In this example, we consider the following network (Figure |4]) where the distance between 
the transmitters (and the receivers) is 8 meters. In this network the interference each receiver 
experiences is lower compared to the received signal. 

The SINR region is convex and hence the required QoS can be easily achieved by simultaneous 
transmissions, even for fixed QoS requirements, as shown in Figure [5| 

The graphs in Figure |6] illustrate the response of our distributed algorithm to the network 
considered. Each column represents one channel. In the first row, the power levels for all users 
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Fig. 4. A wireless ad-hoc network of n — 4 nodes, consisting of two ommunication pairs {Si — >■ Ri}. The distance between 
the transmitters is 8 m. 



SINR Region 




2 4 6 8 10 12 



SINR of user 1 

Fig. 5. The SINR region of wireless ad-hoc network of Figure [4] It is convex and hence any set of SINRs can be realized 
with simultaneous transmissions. 

over time are depicted. In the second row, the allotted SINR targets over time are shown, whereas 
in the third row the figures show actual SINR values for each of the channels. The graphs show 
very fast convergence and since the initial powers are very close, the distributed algorithms 
behave similarly. The system converges to a feasible solution quite fast, even though it is not the 
optimal one in terms of minimizing the individual's total power. In addition, both transmitters 
achieve the required QoS on average. 

B. Example 2: Infeasible network 

In this example, we consider the following network (Figure |7]), which is the same as before, 
but the distance between the pairs is reduced now to 1.8 meters. The interference between the 
two pairs is much bigger now and the SINR region for this network is non-convex (Figure [8]). 
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As shown in Figure [8| the SINR region is non-convex and hence if the SINR QoS requirement 
is outside the SINR region, then each user has to reduce the SINR requirement in the channel 
where the other increases it. In this way, they reduce the interference they cause to each other 
and they can eventually achieve the QoS requirements (Figure |9]). 

VII. Conclusion 

While traditional approaches that tackle the problem of distributed power control consider a 
single channel only, here we investigated the more challenging, yet more promising power control 
problem over multiple channels. More specifically, we have developed a distributed algorithm 
consisting of two update formulae: one for the power and one for the allotted SINR target. In 
this work, we have considered multiple channels and aim to achieve the QoS requirement on 
average over all channels. Using standard Lyapunov stability theory we found the power control 
and allotted SINR target algorithms such that the whole network as a system converges to a 
feasible solution, if one exists. As aforementioned, the solution found is not necessarily the 
optimal. The results of this work are of paramount importance since wireless nodes can make 
use of multiple channels simultaneously and achieve a QoS on average, that otherwise would 
be impossible. Note that heterogeneous networks can be also considered in the sense that each 
node can be equipped with different number of antennae and hence have access to a different 
number of channels. Such algorithms are useful for elastic and/or opportunistic traffic, where 
time-varying rates are allowed and large delays are tolerated. 

The convergence rate of this algorithm is an ongoing research. It is also important to specify 
the stability conditions of this algorithm in the presence of uncertainties (e.g. time-varying delays 
and changing environment). 
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Fig. 6. Power, allotted SINR target and SINR for each communication pair at each channel is shown. The allotted SINR 
converges fast to a feasible solution (not necessarily the optimal) and the system converges quite fast. 
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Fig. 7. A wireless ad-hoc network of n = 4 nodes, consisting of two communication pairs {Si — >■ Ri}. The distance between 
the transmitters is 1.5 m. 



SINR Region 




SINR of user 1 



Fig. 8. The SINR region of wireless ad-hoc network of Figure |7] It is non-convex and if the SINR QoS requirement is outside 
the SINR region, the users have to tend towards different channels. 
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Power Control of individual users in channei 1 



Power Controi of individual users in cfiannei 2 




20 30 
Time (ms) 

SINR of individual users in channel 1 



0.7 
0.6 



0.3 
0.2 
0.1 

C 

5.2 r 

5 - 
4.8 



10 20 30 4 

Time (ms) 

Desired SINR of individual users in cfiannei 2 



4.6 

tr 

-o 4.4 ■ 



° 4.2. 



10 20 30 4 

Time (ms) 

Desired SINR of individual users in channel 1 



20 30 
Time (ms) 



24 
22 
20 
18 
16 

cc 

Z 14 

CO 

12 
10 
8 
6 
4 



20 30 
Time (ms) 

SINR of individual users in channel 2 



20 30 

Time (ms) 



Fig. 9. Power, allotted SINR target and SINR for each communication pair at each channel is shown. The allotted SINR 
converges slowly, since the QoS required cannot be achieved with simultaneous transmission. 
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